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Magneto-optical traps (MOTs) of highly magnetic lanthanides open the door to explorations of 
novel phases of strongly correlated matter such as lattice supersolids and quantum liquid crystals. 
We recently reported the first MOTs of the five high abundance isotopes of the most magnetic atom, 
dysprosium. Described here are details of the experimental technique employed for repumper-free Dy 
MOTs containing up to half a billion atoms. Extensive characterization of the MOTs' properties — 
population, temperature, loading, metastable decay dynamics, trap dynamics — is provided. 

PACS numbers: 37.10.De, 37.10.Gh, 37.10.Vz, 71.10.Ay 
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I. INTRODUCTION 

Ultracold gases of fermionic dipolar atoms provide 
the opportunity to examine — in a pristine and tunable 
setting — the non- Fermi liquid, strongly correlated elec- 
tronic behavior manifest in some of the most interest- 
ing materials of late: high- T c cuprate superconductors, 
strontium ruthenates, 2D electron gases, and iron-based 
superconductors. Strong and competing interactions in 
these systems induce transitions to states beyond the fa- 
miliar insulating, metallic, and superconducting. Specif- 
ically, phases that break rotational and translational 
symmetries emerge in a manner akin to those found 
in classical liquid crystals, e.g., the nematic and smec- 
tic [TJG]. Although quantum liquid crystal (QLC) the- 
ory can describe these non-Fermi liquids in a compelling 
and general framework, unwanted solid state material 
complexity — disorder and dynamical lattice distortions — 
can obscure the underlying electronic physics, and lack 
of wide system tunability can hamper efforts to fully ex- 
plore QLC phase portraits. 

Exploiting large dipole-dipole interactions (DDI) in ul- 
tracold gases will allow the exploration of QLC physics 
in this inherently more tunable and characterizable sys- 
tem. Recent theoretical proposals employing strongly 
magnetic fermionic atoms or polar fermionic molecules 
have begun to shed light on the accessible QLC physics. 
These include predictions of uniaxial (meta-nematic) [3 
and biaxial nematic [4] distortions of the Fermi surface 
in the presence of a polarizing field, and meta-nematic 
and smectic phases in 2D anisotropic optical lattices [5]. 
An exciting prospect lies in the possibility of observ- 
ing spontaneous magnetization in dipolar systems, and 
Refs. [6[ [7] postulate the existence of observable quan- 
tum ferro-nematic phases and spin textures in ultracold 
highly magnetic fermionic atomic gases in zero polarizing 
field. 

While many exciting results will continue to arise from 
the degenerate fermionic dipolar molecule system [5HTT] , 
to most easily observe true (non-meta) QLC phases, the 
symmetries of interest should be spontaneously broken, 
which is not possible when employing ground state polar 
molecules. This is because the strong, r -3 character of 
the DDI is realized only in the presence of a rotational 



symmetry breaking, polarizing electric field that mixes 
opposite parity states. Ultracold chemical reactions can 
also hamper the use of fermionic polar molecules for stud- 
ies of DDI-induced exotic phases in 3D, as recent exper- 
iments in KRb have shown [lOl[TT]. In contrast, highly 
magnetic atoms exhibit the DDI interaction even in the 
absence of a polarizing field and are largely immune to 
chemical and inelastic collisions when spin-polarized in 
an optical dipole trap. While both the electric and mag- 
netic DDI can be continuously tuned to zero, only the 
magnetic DDI can be tuned negative [T2] . 

One must look to the rare earth (lanthanide) series 
to find atoms possessing masses and magnetic moments 
large enough to support exotic bosonic and fermionic 
phases. The extraordinarily large magnetic dipole of dys- 
prosium (10 /i#), which possesses the largest magnetic 
moment of any fermionic atom — and is tied with ter- 
bium in possessing the largest moment among bosonic 
atoms [13] — is likely sufficient to induce QLC phases [5]- 
|7]- We describe here in detail the Dy magneto-optical 
traps (MOTs) [14] that may serve as a precursor to de- 
generate, highly magnetic Fermi — and Bose — gases. 

It is not yet clear whether a degenerate gas of highly 
magnetic lanthanide atoms is possible to create, but over- 
coming a large DDI while cooling to degeneracy has 
precedent in the Bose-Einstein condensation (BEC) of 
chromium [15]. The attainment of Cr BEC opened the 
door to bosonic ultracold dipolar physics [16] . Chromium 
possesses the large magnetic dipole moment of 6 Bohr 
magnetons (//#)• This is 6-times larger than that of the 
alkali atoms' — Rb has a magnetic moment of \i = 1 \±b in 
the doubly-polarized state — and represents a significant 
36-fold enhancement in the DDI strength 
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where r is the distance between two dipoles and 6 is the 
angle between the direction of polarization and the rel- 
ative position between the two particles. However, as 
strong as Cr's DDI may be, recent calculations [17 sug- 
gest that novel lattice phases predicted by the extended 
Bose-Hubbard (eBH) model [16], such as density waves 
and lattice supersolids [18], lie just beyond the reach of 



Cr's capability. Specifically, the DDI energy must dom- 
inate the contact interaction energy, which occurs when 
e = jaoja 2 m/127rh 2 a s > 1, where m is the mass, a s is 
the s-wave scattering length, and the extra factor of 3 
in the denominator designates e > 1 as the regime in 
which homogeneous dipolar BECs become unstable to 
collapse [16]. To observe novel lattice phases, e should 
be >0. 7-0.8 [17]. For 52 Cr, e Cr = 0.15, and even with 
the demonstrated five-fold reduction of a s via a Fesh- 
bach resonance [19 [20 , ec r = 0.7 remains at or below the 
threshold for new phases. In contrast, e^ y = 1.34 which 
is 9x larger than ec r , assuming that the as yet unmea- 
sured scattering length for at least one of Dy's bosonic 
isotopes is approximately equal to 52 Cr's a s = 100 clq. 
With such a large 6D y , exploring supersolids and den- 
sity wave phases without the use of Feshbach resonances 
should be possible with Dy. 



As candidates for fermionic dipolar physics, existing 
MOTs of highly magnetic fermionic 53 Cr and 167 Er are 
not yet populous enough to contemplate cooling to de- 
generacy [21, 22 . The technique of buffer gas cooling has 
been used for cooling lanthanides and has proven success- 
ful in producing large 500 mK samples [23] . While Dy 
has been adiabatically cooled to ~50 mK at final den- 
sities of ~10 9 cm -3 [24 in such an experiment, evapo- 
rative cooling in a magnetic trap to lower temperatures 
has not yet proven effective due to large inelastic colli- 
sions. Thus, the extraordinarily populous 0.1-2 mK Dy 
samples discussed here open the door to a rich landscape 
of physics. 



Ultracold samples of bosonic and fermionic isotopes 
of Dy would find application beyond the QLC and eBH 
physics described above. Co-trapping isotopes of Dy 
present an interesting system to explore Bose-Fermi mix- 
tures of near equal mass, reminiscent of 3 He- 4 He studies, 
but now in the presence of same and cross-species dipo- 
lar interactions. Studies of large-spin degenerate spinor 
gases [25] [26] , simulations of dense nuclear matter [27] , 
creation of unconventional superfluid pairing [28 , and 
explorations of zero sound [29] and roton modes [30] in 
dipolar gases are exciting avenues of research. In ad- 
dition, ultracold samples of Dy will aid precision mea- 
surements of parity nonconservation and variation of 
fundamental constants [31], single- ion implantation [21], 
and quantum information processing [32] [33] . With re- 
gard to the latter, the low-lying telecom (1322 nm) and 
InAs quantum dot amenable (1001 nm) transitions (see 
Fig. HI) will be useful for creating hybrid atom-photonic 
or atom-quantum systems that exploit potentially long- 
lived spin coherences in ensembles of Dy for quantum 
memory. Novel ultracold collisions [34; and complex Fes- 
hbach resonance-induced molecular association phenom- 
ena are expected to appear in traps of these non-6' state 
atoms, which may enable exploration of ID strongly cor- 
related gases with highly magnetic molecules [35] . 
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FIG. 1. (color online). Dy energy level structure [36J - 
138] , where J is the total electronic angular momentum and 
F — J + /, the total electronic plus nuclear angular momen- 
tum. The prime indicates excited state quantum numbers. 
(a) The MOT and Zeeman slower use the strongest laser cool- 
ing transition at 421 nm; there are nearly 140 lower energy 
met ast able states. Several other laser cooling transitions ex- 
ist between the ground state and the excited states. The red 
states are even parity and the odd are shown in black. Dyspro- 
sium has five high abundance isotopes; three bosons ( 164 Dy, 
162 Dy, and 160 Dy with 7 = 0) and two fermions ( 163 Dy and 
161 Dy with I = 5/2). (b-c) The two fermions have oppositely 
signed nuclear spin, resulting in a relative inversion in the 
hyperfine structure. 



II. LASER COOLING SYSTEM 

A simple laser cooling and trapping method for highly 
magnetic atoms has now proven successful for three mag- 
netic lanthanides, Er, Dy, and Tm [14] (2TJ [39] . Measure- 
ments of MOT recycling dynamics show that the working 
principle behind the Er and Dy MOTs are similar [T4]|2T]. 
Despite the existence of more than a hundred levels be- 
tween ground and the open excited state of the cooling 
transition, no repumping laser is necessary because of a 
novel recycling mechanism [21 : After decaying out of the 
MOT cooling transition to metastable states, the highly 
magnetic atoms remain confined in the MOT's magnetic 
quadrupole trap while the metasable state population re- 
cycles to the ground state, at which point the atom are re- 
captured by the MOT. The strongest transition J — » J+l 
transition is used for Zeeman slowing and MOT cooling. 
Whereas weakly magnetic atoms would be lost from the 
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FIG. 2. (color online). Schematic of the Dy laser cooling 
system. See text for details. 
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FIG. 3. (color online). Sketch of the Dy UHV chamber sys- 
tem. The blue arrows depict the slowing, cooling, and imaging 
laser beams. 



trapping region during time in which the atoms spend in 
the metastable states, the strong dipole moments (7 \±b 
for Er) allow confinement even in the MOT's magnetic 
quadrupole gradient. The 1:10 5 branching ratio between 
population decay back to the ground state versus to the 
metastable state is sufficient for a MOT to form given 
the efficient population recycling. 

We highlight four criteria necessary for the large- 
population, highly magnetic, repumperless MOT to work 
when using strong open transitions: 1) to be captured by 
the MOT, the branching ratio must be small enough that 
atoms exit the Zeeman slower in their ground state; 2) 
the decay channel through the metastable states must 
be rapid enough that inelastic collision processes — from 
background, spin-changing, or light induced collisions — 
do not deplete the metastable magnetic trap (MT) reser- 
voir; 3) decay through the metastable MT should be slow 
enough for it to serve as a reservoir to capture more atoms 
than those cycling on the MOT transition; 4) for isotopes 
with / 7^ 0, the difference frequencies between adjacent 
F — >• F' = F + 1 hyperfme transitions in the ground 
and excited states must be small enough t hat th e cooling 
laser itself serves as a repumper (see Sec. IV B). 

All of these criteria are satisfied in the five most abun- 
dant isotopes of Dy, the level structures of which are de- 
picted in Fig. [T] Before discussing the Dy MOT recycling 
dynamics and population in Sec. |IV[ we first describe the 
laser system employed for producing the light necessary 
to drive the Dy laser cooling transition at 421 nm. 

The strongest J — )> J + 1 cycling transition is the 
421-nm line, which has a broad linewidth of 31.9(0.7) 
MHz [36| 140 ] . A total of 2 W of 421-nm optical power 
is generated by two Tiisapphire laser systems. Such 
high power is necessary because the saturation inten- 
sity 58 mW/cm 2 is large and, as we discuss below, the 
Zeeman slower optimally uses 1 W of power. The first 



Tiisapphire system (TiSl) is pumped with 15.2 W at 532 
nm, which when frequency doubled in a ring cavity with 
an LBO crystal produces 600 mW of continuous wave 
(CW) 421-nm light with a 20-kHz linewidth [41]. The 
second doubled Ti:sapphire laser system (TiS2) produces 
up to 1.6 W CW at 421 nm (again using a ring cavity 
frequency doubler) from a 4.3 W 842-nm beam. TiS2 is 
pumped by a 18.5 W 532-nm laser, and the 421-nm light 
has a linewidth of ~50 kHz at 421 nm. 

TiSl produces the 421-nm light for the MOT, trans- 
verse cooling, and imaging beams, while TiS2 is used 
solely for the Zeeman slower. Each laser is locked, at 
842 nm, to its own low-finesse reference cavity. These 
cavities exhibit ~40 MHz/hr frequency drift, so we have 
developed a more stable locking scheme which provides 
a long-term stability of several 100 kHz. The frequency 
reference is derived using the transfer cavity technique: 
the 842-nm light from TiSl is locked to a mode of an op- 
tical cavity whose length is stabilized by a 780-nm laser 
that is resonant with the cavity ~2 x 10 4 modes away. 
This 780-nm laser is itself frequency stabilized to a hy- 
perfme transition in Rb using saturation-absorption spec- 
troscopy. Thus, the stability of a medium- finesse cavity 
locked to a Rb frequency reference is transferred to the 
842-nm laser of the TiSl system. The TiS2 system is 
offset locked [42] at the Zeeman slower laser detuning by 
detecting a beat note at 842-nm between the TiSl and 
TiS2 lasers. See Appendix |A| for more details. 



III. TRAPPING APPARATUS 

Fig. [3] illustrates the schematic of our Dy vacuum 
chamber system. A thermal atomic beam is first gener- 
ated with a high-temperature oven. This atomic source 
is then collimated through a differential pumping tube, 



and a larger flux is achieved by transverse cooling. The 
Dy atoms are decelerated and cooled by a "spin-flip" Zee- 
man slower. Finally, the slowed atoms are captured and 
loaded into a large-gradient MOT in the trapping cham- 
ber. 

The vacuum system consists of four main sec- 
tions: high-temperature oven, transverse cooling, Zee- 
man slower, and MOT trapping chamber. Each section — 
except the Zeeman slower — has an ion gauge to monitor 
the vacuum pressure and a dedicated ion pump. In ad- 
dition, there is a stainless steel tube (inner diameter 4.6 
mm, length 18 cm) that provides differential pumping 
between the oven and transverse cooling sections. The 
Zeeman slower (inner diameter 1.7 cm, length 54 cm) 
also serves as a differential pumping tube between the 
transverse cooling and MOT sections. A titanium subli- 
mator provides additional pumping in the MOT trapping 
chamber. 



A. High-temperature oven and transverse cooling 

Because Dy has a very high melting point of 1412 °C, a 
dual- filament, all tantalum high-temperature effusion cell 
with a water-cooling shroud [43] is used to heat mm-sized 
pieces of Dy. The non-isotopically purified Dy is placed 
in the tantalum crucible and heated to 1250 °C during 
typical MOT operation. A near-uniform temperature is 
achieved by two servos controlling heaters located at the 
the tip and main part of the oven. The crucible with a 5- 
mm diameter orifice generates a Dy beam. By attaching 
a flexible bellow to the oven and controlling it with me- 
chanical positioning stages, the direction of atomic beam 
is optimized for maximum MOT population. 

Because blackbody radiation from the crucible can 
heat the UHV chamber walls, a 2.25" diameter tanta- 
lum shield is installed around the tip of the differential 
pumping tube, immediately before the transverse cooling 
section. Heat propagation through the vacuum system is 
avoided by continuously water cooling the ConFlat (CF) 
flange to which the tantalum disc and differential pump- 
ing tube are welded. The atomic beam may be shut- 
tered by an in-vacuum, pneumatically controlled tanta- 
lum shield positioned a few mm from the crucible orifice. 

Typical oven operation at 1250 °C provides a Dy va- 
por pressure inside the crucible of ~7 x 10 -2 Torr. We 
achieve a vacuum of 1.0 x 10 ~ 9 Torr at the ion gauge of 
the oven section with the use of a 75 L/s ion pump. This 
low pressure is achieved with help from the getter prop- 
erties of the evaporated Dy on the tantalum disk [44] . 
Fifteen grams of Dy lasts ~400 h at 1250 °C operation. 
A gate valve after the oven section maintains the vacuum 
in the rest of system while detaching the oven for the re- 
filling. At the transverse cooling stage, the pressure is 
8.2 x 10 -11 Torr when continuously evacuated with a 55 
L/s ion pump. 

A larger atomic beam flux is achieved by transverse 
cooling the portion of the beam after the differential 
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FIG. 4. (color online). Magnetic field profile of Zeeman slower 
as a function of position away from the center of the trapping 
chamber. Target (line) and measured (o) fields are shown. 
The profile is measured with 3.01 A in Section I, 3.02 A in Sec- 
tion II, and 40.24 A in Section III. Note that Zeeman slower 
was designed for 50% maximum deceleration efficiency. The 
optimal current for MOT production is 1.08 A in Section I, 
2.66 A in Section II, and 40.21 A in Section III. Inset shows 
field zeroing at MOT position. 



pumping tube and before the Zeeman slower. A detailed 
study of Dy transverse cooling appears in Ref. [45 . The 
transverse cooling beams are elliptically shaped such that 
they match the atomic beam; their dimensions are 4.4 
mm by 18 mm. One beam is oriented horizontally to the 
the atomic beam and is retroreflected through the cham- 
ber via ant i- reflect ion coated windows. The second beam 
enters the chamber vertically, passes through a quarter- 
wave plate, into the chamber and is again retroreflected 
through a quarter-wave plate before entering the cham- 
ber. Addition of the quarter- wave plates along the verti- 
cal direction enhances transfer efficiency; no such gain is 
observed for the horizontal branch, presumably since the 
local magnetic field is nonzero. The cooling beams are all 
red-detuned from the 421-nm transition by 0.2-0.4 T. At 
a total power of ~200 mW divided among the horizontal 
and vertical beams, the MOT population is enhanced by 
up to a factor of 4. Additional power slightly reduces the 
MOT population, ostensibly due to increased metastable 
population shelving before th e Zeem an slower (and pos- 
sibly recoil heating). See Sec. IV B for more details. 



B. Zeeman slower 

The spin-flip (zero-crossing) Zeeman slower [46] decel- 
erates and cools the collimated atomic beam from 480 
m/s (most probable velocity) to 30 m/s [47]. The inho- 
mogeneous magnetic field profile of the Zeeman slower is 
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FIG. 5. (color online). MOT population as a function of (a) 
Zeeman slower detuning and (b) Zeeman slower laser power. 
Optimal operation occurs for a detuning of -21 T. Population 
saturates above 1 W of laser power. 



shown in Fig. [4j and the positive and the negative fields 
are generated from three distinct coils. The first section, 
which is closest to the oven, has 11 staggered layers of 
magnet wire helically wound around the 1.1" CF nipple 
forming the Zeeman slower. The second section has 18 
layers with current flowing in the opposite direction to 
the first. A third section is added to cancel the stray 
field from the first two sections at the MOT position. It 
consists of a 10- turn coil of hollow square magnet wire 
(4.2 mm on a side) through which chilled water circulates. 
Because of the heat dissipated in the first two solenoid 
sections, chilled water through copper tubing wrapped 
around them is necessary to keep the temperature below 
30 °C. 

Up to 1.5 W of Zeeman slower beam power (Fig. [5| is 
used to slow the atomic beam. The laser is detuned -21 T 
[Fig.p)Fa)] from the 421-nm atomic transition [48] and the 
beam is focused into the orifice of the oven crucible (see 
Fig. [3|. A 45° aluminum mirror is installed inside the 
vacuum system to avoid Dy coating the entrance view- 
port of the Zeeman slower laser beam. The Dy film coat- 
ing the mirror is reflective at 421 nm. The diameter of the 
slowing laser is 2.5 cm at the vacuum window entrance. 
As shown in Fig. |5|b), total trapped atom population 
^Vxotai saturates at ^1 W of input power; presumably 
power broadening aids the velocity capture range of the 
slower. 



vide access for the atomic beam and counterpropagating 
Zeeman slower beam. 

The magnetic quadrupole field used for both the MOT 
and magnetic trap [14] is generated from a coil pair 
in near-anti-Helmholtz configuration. The coils gener- 
ate V 2 B = 0.69G/cm-A along i, which points along 
the quadrupole axis of symmetry; each coil has a cross- 
section of 10 rows (z) and 7 columns (p). The coils are 
water cooled to support the ~30 A used for the MOT. 
The electric current is controlled by a servo providing a 
440 {is turn-off time. Stray field cancellation coils reduce 
the residual field to < 1.8 mG. 

We achieve a MOT chamber pressure of typically 
1.2 x 10 -11 Torr during MOT operation. The vacuum 
is created with the help of a 75 L/s ion pump and a ti- 
tanium sublimation pump. We have noticed a lowering 
of the vacuum pressure from the evaporated Dy after the 
oven has run for several hours. 



D. Imaging 

After the atoms are released from the trap, absorption 
images are taken from the CCD camera installed on the 
side of the trapping chamber with a probe beam pulse 
width of 60 /is. This short pulse width ensures that 
that atoms are approximately stationary and shelving 
into metastable states does not occur during the imaging 
process. In order to take reliable temperature measure- 
ments of the MOT while avoiding eddy currents (lasting 
~1.5 ms) incurred during magnetic quadrupole coil shut- 
off (at td — 0), the time-of- flight data were taken in 500 
/is steps from td — 3 to 5-7 ms. 



IV. DY MOT CHARACTERISTICS 

In this section we discuss a detailed, semiclassical 
model and supporting data for how the repumper-less 
Dy MOT functions. Measurements of the Dy MOT pop- 
ulation and temperature versus MOT loading time, in- 
tensity, detuning, and magnetic quadrupole field gradient 
are presented. We also compare Dy MOT characteristics 
to those of the Er MOT [2X1151]. 



C. MOT trapping chamber 

The MOT employs a three-retroreflected-beam config- 
uration (a six-beam configuration can also be used, with 
similar results). Each beam is aligned and collimated in 
free-space with a waist [49] of 1.1 cm, and the MOT typ- 
ically has a detuning of A = —1.2 T. A total intensity 
of 0.17-0.2/ s provides the maximum MOT population, 
where I s ~ 2.7 x 58 mW/cm 2 [50 . A stainless steel oc- 
tagon chamber with two 6" diameter CF viewports on 
top and bottom and six 2.75" CF viewports on the side 
provide the optical access necessary for the MOT beams 
and imaging beam. The remaining two CF ports pro- 



A. MOT decay and motional dynamics 

As mentioned in Sec. [TT] the Dy MOT functions with- 
out any repumpers due to a population recycling mech- 
anism based on its extraordinarily large magnetic mo- 
ment. Figure [6] depicts a semiclassical rate equation 
model describing this recycling process. Population is 
lost from the excited state of the MOT at rate f e ^R\ — 
which increases for larger excited state fractions / ex — and 
decays to one of the (possibly many) metastable states 
(see Fig. [I]). A fraction of atoms 1 — p are either in the 
m = Zeeman substate or in a strong magnetic field 
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FIG. 6. (color online). Dy MOT recycling and continuously 
loaded MT schematic. 



seeking m state and never become confined in the MT. 
The remaining atoms decay either through a fast decay 
channel at rate i?f as t and fractional population 1 — q or 
through a slow channel at rate R s \ w Atoms whose elec- 
tronic population reaches the ground state are reloaded 
into the MOT at rate i? r eioad> which depends on the total 
MOT beams' saturation parameter, s = 7/(l + (2A/T) 2 ), 
where I = I/I s and A is the detuning from resonance. 
There must be some loss R\ OSS 2 from the magnetic trap as 
the population cascades, but as shown in Ref. [14 , spin 
relaxation loss is 20 x slower than the timescales of the 
measured i?f as t and R s \ow 

The diagram in Fig. [6] may be represented with the 
following rate equations: 

A^MOT = ^reload ^MT ~ Jex-Rl^MOT + i?L> 
^Vfast = (1 - q)pfe*RlN M OT ~ ^fast^fast, 
^slow = QPfe^Rl^MOT ~ ^slow^slow, 
7V M T = ^fast^fast + ^slow^slow ~ ^reload^MT, (2) 

where A/"mot, ^Vf ast , N s \ ovn and Nmt are the populations 
of the MOT, fast (slow) metastable state decay channel, 
and ground state MT, respectively, and R\, is the loading 
rate from the Zeeman slower. Fits of these rate equations 
to MOT population decay curves allow the determination 
of all free parameters for both the 164 Dy boson and the 
163 Dy fermion: 



[R\ 
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pl63 pl63 1 
^fast > ^slowj 



[1170(20), 19(2), 1.5(1)] s~\ (3) 



[^r,«^si 6 o 4 w] = [1700(100), 29(1), 2.3(1)] s" 1 , 
\p,q] = [0.82(1), 0.73(1)], 

where p and q for the two isotopes are equal within error. 
An example of the Dy MOT decay and fit are shown 
in Fig. [7^a). To obtain these data, we collect MOT 
fluorescence using a 200-mm 2" aspherical lens and fo- 
cus it through a pinhole, using a 60-mm 2" aspherical 
lens. The light is then collected on an amplified, high- 
bandwidth PIN photodetector. We load the MOT for 
2-3 s to reach steady-state population before extinguish- 
ing, at t m = 0, the MOT beams, the Zeeman slower 
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FIG. 7. (color online), (a) Population ratio A^mot of re- 
captured MOT to steady state MOT. Black line is fit using 
Eqs.|2]to iVMOT with Zeeman slower and atomic beam off, and 
t m = Is delay between steady state MOT and recapture, (b)- 
(c) Oscillations in fluorescence at the peak of the recaptured 
MOT population appear for MOT beam saturation parame- 
ters s< 1. (d) The Dy MOT operates in the mechanically 
underdamped regime, as shown by breathing mode oscilla- 
tions, and the photon scattering rate seems to be correlated 
with cloud diameter. 



beam, and the atomic beam (via the in- vacuum shutter) . 
Note, the magnetic quadrupole field remains on through- 
out this experiment, providing confinement for weak- field 
seeking metastable and ground state atoms. We wait for 
t m = 1 s before recapturing the MOT by turning on 
the MOT beams, but we leave the Zeeman slower beam 
and atom beam off to ensure there is no external load- 
ing of the MOT during the measurement, i.e., Rl = 0. 
Shorter t m times bias the results due to ill-defined initial 
conditions: we fit the data assuming nearly all popu- 
lation has decayed to the ground state, allowing us to 
set N MOT = 7V fast = N s}ow = at t - t m = 0, while 
TVmt is left free to vary and is normalized to the steady 
state MOT fluorescence signal. Fluorescence traces as in 
Fig. ma) are averaged over 16 sequential runs on a digital 
oscilloscope. 

We noticed an oscillation of the recaptured MOT fluo- 
rescence near peak signal when operating at small satura- 
tion parameter s«l. The pinhole does not restrict nu- 
merical aperture so severely that the modulation is sim- 
ply due to a loss of photons when the cloud is large. This 
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FIG. 8. (color online). Fits of Eq.|5]to f ex Ri versus I = I/I s 
and A for both fermionic 163 Dy and bosonic 164 Dy. A simul- 
taneous fit to the nine (eight) data points for 163 Dy ( 164 Dy) 
allow the extraction of Ri from the product f ex Ri via Eq. [5] 



oscillation disappears when the intensity of the MOT 
light is increased and/or the detuning of the beams de- 
creased to the point that s ~ 1 [compare Figs, fflb) and 
(c)]. Data corrupted with oscillations could not be fit 
using Eqs. [2j so all the data used to extract R±, i?f ast , 
^siow. P, and q were taken near s ~ 1. 

A possible explanation of these fluorescence oscilla- 
tions lies in the motional dynamics of the Dy MOT. Most 
MOTs, e.g., of Rb, are solidly in the overdamped regime, 
meaning they have a ratio a — Tmot/^mot > 1 of the 
MOT's optical damping rate Tmot = 13/ 2m to the trap 
oscillation frequency ujmot = y/ft/m. In a ID treatment, 
the damping coefficient /3 is 



/3 



8M 2 |A|7 



r(i + / + (2A/r) 2 ) 



2^2' 



(4) 



where k is the wavenumber of the MOT light, mass is 
m, and the spring constant is n = n'VBfi /hk [46]. In 
this latter expression, \jl = (g e m e — g g m g )iiBi where gi 
and rrii are the g factors [52] and Zeeman substates of 
the ground and excited levels. For a typical Rb MOT 
operated to maximize trap population, a^b ^ 7. For the 
Er MOT with the parameters used in Ref. [21], «Er = 1.3 
and oscillations where not reported [53]. However, the 
Dy MOT, when operated in the regime that maximizes 
trapped atom population (see Sec. IV E), is just in the 
underdamped regime ar> y = 0.8. The origin of this lower 
a lies in the combination of the larger (smaller) m, &, 
T, // and VB (7 and A/r) parameters in the highly 
magnetic Dy MOT versus typical Rb MOTs. 

Indeed, a& y = 0.8 corresponds to an oscillation pe- 
riod of 4 ms, and a measurement of the breathing mode 
of the Dy MOT, shown in Fig. [7Fd), shows a damped 
oscillation with a period of ~6 ms after the magnetic 
quadrupole field gradient is decreased by a factor of 2 at 
t m = 0. We find that the breathing mode period scales 
oc yVB. The similar periods support the notion that 
fluorescence is modulated by the motion of the atoms 
during recapture. (Atoms confined in the magnetic trap, 
from which the MOT is recaptured, possess a different 
spatial profile from those in the MOT.) The breathing 



mode should have a period half that of the trap oscil- 
lation frequency, indicating that the actual a^ y is just 
slightly below unity; a 3D numerical MOT calculation in- 
corporating magnetostatic forces on the highly magnetic 
Dy could better estimate a for these system parameters. 
We conjecture that optical pumping, Zeeman shifts, and 
spin polarization near the cloud edge — as well as modula- 
tion of optical density — could change the effective photon 
scattering rate as the cloud expands and contracts. To 
avoid complications caused by these motional dynamics, 
we operate the MOT in the overdamped regime of s ~ 1 
for the decay rate measurements. 

To extract MOT decay rates R\ from the product 
fexRi — as wen as to obtain better error estimates of the 
other model parameters — we repeat the measurement in 
Fig. ma) for several combinations of 7 and A. MOTs of 
bosonic 164 Dy and fermionic 163 Dy are studied in a sim- 
ilar manner to investigate differences in decay dynamics 
due to hyperfine structure. Figure [8] shows the set of data 
simultaneously fit to the function 



R(I, A) = R 1 I/(2 + 21 + 2(2A/r) 2 



(5) 



where R = f ex Ri- Using this method to determine R\, 
we arrive at the rates in Eq. Q for these two Dy isotopes. 

Detailed simulations of the decay channels still seem 
beyond the reach of tractable calculations, but such cal- 
culations might be able to employ the following simplifi- 
cation for modeling the fast decay channel. The shortest 
possible decay channel (the fast channel), must involve 
at minimum two metastable levels due to the need to 
switch parity from odd (421-nm exited state) to even 
to odd before decaying back to the even parity ground 
state. Recent calculations reported in Ref. [54] indicate 
that the lifetimes for the 1001-nm and 1322-nm levels 
are 2i\ • [53,23] s _1 , respectively, which indicates that 
rapid decay to these levels followed by a delay given by 
their lifetime is a likely candidate for the 20-30 s _1 fast 
decay channel. For the slow decay channel, it is pos- 
sible that the population becomes shelved in a longer- 
lived state that decays to ground or an intermediate state 
via a non-electric dipole (El) allowed transition. Small- 
energy-difference decays could also contribute to the slow 
channel. 

The large q parameter — 73% of the atoms decay 
through the slow versus fast channel — serves a very use- 
ful purpose in providing a continuously loaded magnetic 
reservoir of metastable atoms for the MOT. Indeed, data 
presented in Ref. [14 show that 2.5 x ( 164 Dy) to 3.5 x 
( 163 Dy) more atoms are in the metastable MT plus MOT 
than in the visible MOT alone, resulting in a maximum 
number of laser cooled and trapped atoms nearing 5x 10 8 
for the Dy system. This is reminiscent of the continuously 
loaded MT in the Cr system [55] that provides sufficiently 
large samples, 10 8 atoms, to enable Bose-condensation of 
Cr US]. 
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FIG. 9. (color online). Comparison of relative MOT popula- 
tion to natural isotope abundance for the five highest abun- 
dance isotopes of Dy. Relative population is computed as the 
fraction of an isotope's MOT population out of the sum total 
MOT population for all the isotopes. The fermions, 163 Dy 
and 161 Dy are highlighted in grey. The green diamonds are 
the MOT populations for the fermions when multiplied by 
6/5 for 163 Dy and 6 for 161 Dy to account for poor optical 
pumping of the six hyperfine states during laser cooling. See 
text for details. 



B. MOT population versus isotope 

Figure [9] plots the maximum MOT population ob- 
tained for each stable isotope (except for low 0.1% abun- 
dance 158 Dy) along with the natural abundance of the 
isotope. To facilitate comparisons to natural abundance, 
each MOT population is normalized by the sum of all 
the MOT populations of the isotopes (for the fermionic 
isotopes, the MOT population is adjusted for hyperfine 
shelving; see below). The bosonic MOTs have relative 
populations very close to their proportion of the natural 
abundance. However, the fermionic isotopes 163 Dy and 
161 Dy have populations that are 84.0% and 17.3% their 
natural abundances, respectively. We try, in the follow- 
ing, to argue why repumper-less fermionic MOTs form 
and are observed with these population ratios. 

Unlike the bosons, the two fermions have non-zero nu- 
clear spin / = 5/2, and the opposite sign of the nuclear 
magnetic moment between the two isotopes results in op- 
positely ordered hyperfine levels with respect to F versus 
energy [see Fig. [I] (b) and (c)]. With a total electronic 
angular momentum of J = 8, the isotopes have six hy- 
perfine states, F = 11/2 to 21/2 and F' = 13/2 to 23/2, 
in their ground and 421-nm excited states, respectively. 
Without hyperfine repumpers (a repumper is necessary 
for Rb and Cs MOTs), one would expect that no MOT 
could form on the F = 21/2 —> F' = 23/2 cycling transi- 
tion due to rapid decay to F < 21/2 states. Alkalis such 
as Rb and Cs, however, have smaller hyperfine splittings 
(<270 MHz) between their highest F' excited states than 
does 163 Dy, whose splitting is 2.11 GHz. However, such 
an explanation — that scattering to the lower F hyper- 
fine state is slower due to the larger detuning — ignores 
the 6x shorter lifetime of the Dy excited state. Together 
this implies a depumping rate similar to Rb and Cs. The 
crucial difference, however, lies in the fact that the hy- 



perfine splittings in the ground and excited states of Dy 
are nearly matched, whereas for Rb and Cs, they are 
different by more than 25 x. Thus for Dy, the trapping 
light near-resonant with the F = 21/2 -)► F' = 23/2 cy- 
cling transition also serves as an efficient repumper for 
the F — 19/2 — >> F r = 21/2 transition, thereby prevent- 
ing complete population shelving into dark states. 

A full modeling of the intra- hyperfine manifold popula- 
tion transfer as the atoms progress through the transverse 
cooling, Zeeman slower, and MOT stages is beyond the 
scope of this paper, but we can make a simplified esti- 
mate of the F = 21/2 state decay and repumping rates in 
each of the stages. The atoms spend roughly 20 jjls in the 
transverse cooling beams, scatter ~8 x 10 4 photons from 
the F' — 23/2 state in the Zeeman slower, and spend, 
on average, 25 ms in the MOT scattering photons from 
the 421-nm excited state before decaying to metastable 
states. Accounting for laser intensities (power broaden- 
ing) and detunings, this implies that for 163 Dy atoms, <1 
photon is off- resonantly scattered on the F' = 21/2 state 
(which can decay to F = 19/2) in the transverse cooling 
stage, while nearly 100 and 50 are scattered in the Zee- 
man slower and MOT stages, respectively. However, the 
repumping rates in the Zeeman slower and MOT stages 
back to the F' = 21/2 state from F = 19/2 are only a 
factor of two smaller than the cycling transition scatter- 
ing rate and are 570 x and 1300 x faster than the decay 
rates to F = 19/2, respectively: in very little time the 
population is repumped back into the cycling transition. 
No secondary laser is necessary to repump the 163 Dy sys- 
tem. 

The situation is different in the 161 Dy system due to 
the hyperfine structure inversion and the nearer detun- 
ing of hyperfine levels at the larger F and F' end of the 
spectrum. In the transverse cooling, Zeeman slowing, 
and MOT stages, approximately 10, 3 x 10 3 , and 2 x 10 3 
photons are off-resonantly scattered to the F r = 21/2 
state, respectively; the 161 Dy system off-resonantly scat- 
ters more than in the 163 Dy system. Nevertheless, the 
repumping rates are within a factor of two of the depump- 
ing decay rates, and atoms are repumped nearly as fast as 
they are depumped from the cycling transition. Again, 
a repumper-less MOT is able to be formed with 161 Dy 
despite its multitude of hyperfine levels. Reference [21] 
reported a fermionic 167 Er MOT (I = 7/2), which we 
suspect works in a similar manner. 

This leaves open the question why we observe as much 
as 84% of the natural abundance of 163 Dy, while as little 
as 17% of 161 Dy. The hyperfine population in Dy arriving 
from the high-temperature oven is likely to be distributed 
with a bias toward the higher- multiplicity large-F states. 
This may contribute to the large MOT population in 
163 Dy, but 161 Dy's much lower population may be due 
to the competing effect of less-efficient optical pumping 
and repumping to the F = 21/2 state. 

Optical pumping of population into the F = 21/2 state 
must occur before the Zeeman slower detunes the non- 
F = 21/2 state atoms away from resonance. Indeed, we 
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FIG. 10. (color online). Typical MOT loading data as a func- 
tion of time (blue, solid curve). The trapped atom population 
is detected by a photodetector recording MOT fluorescence. 
The Zeeman slower laser and atom beam shutters open at 
t — 0. The green (long-dash) curve is a simulation using 
Eqs.[2]and[3] 



have already seen that the 421-nm laser tuned to the 
F = 21/2 — »> F' = 23/2 transition can pump atoms out 
of the F = 19/2 state on relatively short time scales, and 
we now estimate all the relative F — >> F' = F + l effective 
scattering rates to examine the repumping efficiency from 
the F < 19/2 states: 
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These ratios are proportional to the number of photons 
scattered from the F = 21/2 — >• F' = 23/2 cycling transi- 
tion, denoted F 91 ,1 and normalized to F}?% = F, 



21/2 



21/2 



161 
'21/2 



1000, during the time the atoms transit the transverse 
cooling stage and first 5 cm of the Zeeman slower [56] . 
F ii 6 /2 = 3 for the F = U / 2 -^ F ' = 13 / 2 transition cor- 
responds to roughly 20 scattered photons in our experi- 
ment, assuming that the population is already optically 
pumped to the m g = F state [57]. This is merely an 
upper bound to the number of photons a non-polarized 
sample would scatter since m g < F — >• m e < F + 1 tran- 
sitions have smaller transition strengths than the cycling 
transition m g = F —> m e = F + 1. 

A striking difference between the isotopes is evident. 
The single 421-nm F = 21/2 -)► F' = 23/2 laser scatters 
many more photons from F < 21/2 states in 163 Dy than 
in 161 Dy. This leads to much more efficient optical pump- 
ing to the F = 21/2 state in 163 Dy. We conjecture that 
this is why 84% of the 163 Dy are trapped but only 17% of 
the 161 Dy. Coincidentally, these percentages are close — 
within experimental error — to 5/6 and 1/6, respectively. 
It is tempting — but probably an oversimplification — to 
conclude that for 163 Dy, all levels but the lowest F state is 
optically pumped to F = 21/2. The case of 161 Dy is cer- 
tainly more complicated: while it is likely that F < 17/2 



states are never efficiently pumped to F = 19/2 (let alone 
F = 21/2), the initially larger F population from the 
oven should lead to larger populations than 17%. A pos- 
sible explanation lies in the inversion of the hyperfine 
state energy hierarchy between the two isotopes. This 
causes Doppler shifts in atoms not optimally decelerated 
in the Zeeman slower to move closer to (further from) 
resonance with the F < 21/2 levels in 163 Dy ( 161 Dy). 
This contributes to a more efficient slowing — and there- 
fore more efficient MOT loading — of 163 Dy versus 161 Dy 
atoms. 

A thorough optical pumping simulation including all 
the hyperfine levels and the actual experimental param- 
eters and geometry could better address these questions. 
A repumping laser system is under construction to iden- 
tify into which F < 21/2 states population accumu- 
late, and we suggest that only one repumper on the 
F = 11/2 H> F' = 13/2 transition for 163 Dy would be 
necessary to capture the full fraction of this isotope's 
natural abundance in the MOT. Reference [45] presents 
other repumping schemes. 



C. MOT loading 



Figure [10] shows a typical MOT population loading 
curve after the Zeeman slower laser and atomic beam 
are unblocked at t = 0. The 164 Dy MOT population- 
proportional to the MOT florescence (atoms in the 
metastable MT remain dark) — rises to a steady state of 
2.5 x 10 8 within 50 ms (blue line). A simulation of Eqs.[2] 
predicts, however, that within 5 s the population should 
reach a steady state 5x larger (green dashed line). The 
Dy MOT population limit could arise from light-induced 
two-body collisions as in the Cr MOT [59 . Suppression 
of the i^reioad term in Eqs. [2] mimics the experimental 
data for a loading rate of R\, ~ 10 10 s _1 . Reproduction 
of the loading data requires a small i? re ioad m simulations, 
while fits to decay data such as in Fig.[7Ja) require i? re ioad 
to be of the order 10 3 s _1 . A main difference between the 
MOT loading and MOT decay experiments is the pres- 
ence of the 1 W Zeeman slowing beam, which could en- 
hance light-induced losses of Dy and reduce i? re ioad- An- 
other possibility involves 2-photon ionization loss while 
the atoms are in the metastable states (twice the cooling 
light energy is close to the ionization potential) [58] . 



D. Comparison to Er MOT decay dynamics 

As discussed earlier, the Er MOT did not exhibit un- 
derdamped oscillations in MOT decay florescence, which 
we believe is due to its smaller magnetic moment. Other 
differences between the Dy and Er MOT decay dynamics 
are presented in this section. 

The analysis of the Er MOT decay dynamics [21] ac- 
counted for loss only via an i^i OS sMT term (i.e., p = 1) 
and with only one decay channel through the metastable 
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states (q = 0) [60 J. Allowing q to be nonzero (thus intro- 
ducing two decay channels), provides a much better fit to 
the Dy MOT decay data. The Er apparatus had a 1000 x 
worse vacuum pressure, which might have obscured the 
long-time tail of the MOT decay and prevented a mea- 
surement of a second metastable state decay channel. 

Perhaps coincidentally, the bosons 164 Dy and 168 Er 
have the same, within experimental error, rate of decay 
to metastable states ~1700 s _1 and branching ratios. In- 
terestingly, the fermion 163 Dy, which possesses hyperfme 
structure, has a significantly smaller decay rate 1170 s _1 
as well as shorter metastable decay rates compared to 
164 Dy. While it is not clear why all the 164 Dy decay rates 
are ~1.5x larger than those for 163 Dy, the existence of 
hyperfme structure in 163 Dy likely plays a role. 

As shown in Ref. [14 , 163 Dy's smaller decay rates 
translate into larger steady state metastable MT pop- 
ulations, which help to enhance its total trapped popu- 
lation. The large q and small R s \ow that provide popula- 
tion enhancement in the metastable MT was also seen in 
MOT-loaded magnetic traps of 168 Er, albeit to a lesser 
extent due to poor lifetimes from a high vacuum [51] . The 
majority of Dy atoms decay through a slower metastable 
channel (2.3 (1.5) s _1 versus 4.5 s _1 for 168 Er). However, 
27% decay through the fast channel at rates 6x (4x) 
larger than R s \ ow in 164 Dy ( 163 Dy). From simulations of 
Eqs. 1 in Ref. [21 and Eqs. [2] here, one can see that the 
measured rates imply that the Er metastable MT reser- 
voir could hold several times more atoms as Dy's; how- 
ever, the severity of Er's magnetic trap inelastic loss rates 
are not yet known. Dy collisions limited the metastable 
trap population to 80% [14] of the maximum achievable 
given by the rates in Eqs. [3] 
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FIG. 11. (color online), (a) Time-of- flight absorption image 
of 164 Dy MOT at t d = 3 ms; MOT fields are extinguished at 
td = 0. Intensity integrations along (b) z and (c) p-directions. 
Grey arrow points to region of minority component for this 
MOT. (b) Double-Gaussian fit (white line) to MOT expan- 
sion (the red line is a single- Gaussian which results in a poor 
fit). The majority component (hot outer cloud) is denned as 
the atoms contributing to the broader Gaussian, while atoms 
in the narrower Gaussian comprise the minority component. 
(c) The anisotropically cooled stripe hampers majority com- 
ponent temperature measurements in the z-direction. The 
origin and temperature characteristics of the minority com- 
ponent are studied in Ref. [51] . 



E. MOT population and temperature 

As reported in Ref. [14], when the power between the 
three sets of MOT beams are mismatched, the Dy MOT 
can exhibit a sub-Doppler cooled core, which depending 
on the details of the power balancing, can be anisotropic 
in its temperature distribution. However, the number of 
atoms in this ultracold core is typically less than 10% of 
the total MOT population. We leave the study of this 
minority component to Ref. [61] , and instead present here 
the temperature and population of the majority compo- 
nent of the Dy MOT as a function of MOT beam inten- 
sity J, detuning A/T, and quadrupole gradient VB Z (see 
Fig. 11 ). Data presented here are typical of the majority 



component in all Dy MOTs. 

To distinguish the hotter majority from the colder 
minority component — when it exists — we fit a double- 
temperature distribution to the time- of- flight expansions 
of the MOT (Fig.[llpb) and Fig. 4 in Ref. [14]). Figure[l2 
shows a compilation of temperature and population data 



of the majority component. In both Figs. 11 and [12] the 
data were taken for a ratio of MOT powers in the z versus 
p directions of I z /I p « 1 and I > 0.2 (J = I Z + 2I P ). This 



results in a MOT with a vertically oriented, anisotropi- 
cally sub-Doppler cooled core. We choose to discuss the 
majority component data with the MOT beams set to 
Iz/Ip < 1 an d I > 0.2 since it is in this regime that we 
obtain the most populous MOTs. 

The major component temperature T ma j or i ty data in 
the p-direction are shown along with plots of Eq. [7| 
the temperature Td from simplified Doppler cooling the- 
ory [46]: 



nv i + / + (2A/r) 2 



4fe B 



2|A|/r 



(7) 



The Doppler cooling limit for Dy on the 421-nm transi- 
tion is 770 /iK, and for the I and A/T's typically used in 
the experiment, the Doppler cooling temperature is ~1 
mK. The prolate MOT T ma j or i ty data are close to that 
predicted by Doppler cooling theory for large s (large 7, 
small A/r) and high gradients. However, T ma j or i ty of the 
atoms at small s and small gradient is below the Doppler 
limit. 

The low gradient and low saturation data are consis- 
tent with the several hundred /iK temperatures of the 
Er MOT [51 . One-dimensional numerical sub-Doppler 
cooling simulations presented in Ref. [51 indicate that 
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FIG. 12. (color online), (a-c) Temperature characterization 
of the majority component of the Dy MOT versus intensity 
I — I /I s , detuning A/r, and MOT magnetic field gradient 
V 2 B (V P B w V z B/2). Data are the temperatures of the 
MOT in the p direction (quadrupole plane of symmetry, grav- 
ity points in the — z direction). Light blue curves are plots 
of Eq. u\ for the following parameters: (a,d) A/r = —1.2, 
V^B = 20 G/cm; (b,e) I = 0.17, V^B = 20 G/cm; (c,f) 
I = 0.20, A/r = -1. (d-f) Total number of atoms (major 
plus minor component) versus /, A/r, and V Z B. The typical 
visible MOT maximum population is A^ max = 2 x 10 8 . Max- 
imum density of ^10 10 cm -3 occurs for: / = 0.33 in (a,d); 
A/r = -1 in (b,e); and V Z B = 20 G/cm in (c,f). 



population-wide MOT sub-Doppler cooling arises from 
the near-equal Lande g factors in Er's ground and ex- 
cited states A#Er = 0.004, which means that the MOT 
magnetic field Zeeman shifts the ground and excited mj 
levels by nearly equal energies. Spherically symmet- 
ric intra-MOT sub-Doppler cooling has been observed 
in the cores of MOTs of more commonly used atoms 
(e.g., Rb [62 J, but unlike Er, the majority of atoms are 
at hotter Doppler temperatures. The Dy MOT — which 
also possesses a near-degeneracy of g factors on the 421- 
nm line [63] — also exhibits population-wide MOT sub- 
Doppler cooling like the Er MOT, which is confirmed in 
ID numerical simulations p~4[[53|. Unlike Er, however, 
the Dy MOT population assumes a double-temperature 



distribution at large s and V Z B. 

In contrast to the Er MOT, we have also ob- 
served anisotropic sub-Doppler cooling to much lower 
temperatures — down to 10 jaK — in the central cores 
of the Dy MOTs. This phenomenon of ultracold 
stripes, shown in Fig. [TT] disappears at low MOT 
beam intensities — coincident with the appearance of sub- 
Doppler cooling in the majority component of the Dy 
MOT. The minority stripe components are typically 1- 
10% the population of the entire MOT, and as such are no 
more populous than 10 6 to 10 7 atoms (see Fig. 11). The 
Er MOT presented in Ref. [51] contained no more than 
2 x 10 5 atoms, and the largest Er MOT to date contained 
1.6 x 10 6 atoms [2T] , We conjecture that the Er popula- 
tion was limited by low Zeeman slower power, which was 
<100 mW [53]. More populous Er MOTs might also ex- 
hibit double temperature distributions when larger MOT 
beam power is employed. The origin of this anisotropi- 
cally sub-Doppler cooled minority component located in 
the Dy MOT core is explored in detail in Ref. [61]. 



V. CONCLUSION 

We presented a detailed description of the creation of 
large-population open-shell lanthanide MOTs with the 
repumper-less technique. The demonstration that the 
repumper-less MOT technique works for Dy in a manner 
similar to Er (but with important differences) lends sup- 
port to endeavors aimed at forming MOTs of other inter- 
esting highly magnetic lanthanides such as holmium [33] . 
We also presented an in-depth characterization of the 
Dy MOT population, temperature, loading rate, isotope 
trapping efficiency, mechanical dynamics, and metastable 
state recycling dynamics. Future work will attempt to 
narrow-line laser cool Dy on the 741-nm transition in a 
manner similar to that demonstrated for Er [64] . 

We note that Dy MOTs made in a similar manner 
to alkaline-earth MOTs [65, 66 — i.e., using the 421-nm 
light for Zeeman slowing (and transverse cooling) but 
creating the MOT on the closed 140 kHz-wide 598-nm 
transition [36] — will likely be possible and would pro- 
vide lower initial temperatures. However, for maximizing 
phase-space density, it is not yet clear whether a lower 
initial temperature outweighs the loss of the metastable 
MT reservoir and associated population gains. It remains 
to be seen whether such MOTs could also trap fermionic 
isotopes without additional repumpers. 
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Appendix A: Transfer cavity lock 

We provide here further details on the transfer cavity 
lock, which is depicted in Fig. [2] The external cavity 
diode laser at 780 nm is locked to a Fabry-Perot cavity 
using the Pound-Drever-Hall (PDH) technique [67 j. The 
Fabry-Perot cavity consists of piano and concave (r = 25 
cm) mirrors with high reflectivity (99.9%) at both 780- 
nm and 842-nm, and the mirrors are attached to a fused 
silica spacer with a length of 10 cm. In order to scan 
the resonance frequencies of the cavity, a piezoelectric 
transducer (PZT) ring is placed between one side of the 
spacer and the concave mirror. With the measured cav- 
ity length, the corresponding free spectral range is 1.4 
GHz, and the cavity has a linewidth of 1.4 MHz. The 
cavity is isolated from acoustic vibration and placed in- 
side a vacuum chamber. The laser is locked to the cavity, 
which provides short time scale frequency stability, and 
then the cavity length is adjusted to bring the cavity and 
laser in resonance with a hyperfine transition of 87 Rb. 
Feeding back to the cavity PZT an error signal derived 



by Rb saturation absorption spectroscopy ensures that 
the cavity resonance frequencies drift by no more than a 
few hundred kHz. 

The cavity now provides a stable frequency reference 
for TiSl's 842 nm beam, which is picked-off from the 
main beam that is directed into the ring cavity fre- 
quency doubler. To bridge the frequency gap between 
the nearest 842 nm mode and the doubled 842-nm wave- 
length necessary for generating the MOT light, the 842- 
nm beam passes through an optical fiber-based electro- 
optical modulator, which imprints strong sidebands at 
a frequency adjustable between 1.4 GHz-2.8 GHz. A 
phase-locked loop generates this microwave frequency 
from a stabilized and tunable RF source. One of the 
sidebands is adjusted in frequency to be resonant with 
the optical cavity, allowing the 842-nm to be locked to 
the cavity using the PDH technique. Once locked, the 
421-nm light from both TiSl and TiS2 (which is offset 
beat note locked to TiSl) may be scanned to the correct 
frequency by adjusting the RF reference. 
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